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ABSTRACT 

In several gamma-ray bursts (GRBs) excess emission, in addition to the standard 
synchrotron afterglow spectrum, has been discovered in the early time X-ray observa- 
tions. It has been proposed that this excess comes from black body emission, which 
may be related to the shock break-out of a supernova in the GRBs progenitor star. 
This hypothesis is supported by the discovery of excess emission in several GRBs with 
an associated supernova. Using mock spectra we show that it is only likely to de- 
tect such a component, similar to the one proposed in GRB 101219B, at low redshift 
and in low absorption environments. We also perform a systematic search for black 
body components in all the GRBs observed with the Swift satellite and find six bursts 
(GRB 061021, 061110A, 081109, 090814A, 100621A and 110715A) with possible black 
body components. Under the assumption that their excess emission is due to a black 
body component we present radii, temperatures and luminosities of the emitting com- 
ponents. We also show that detection of black body components only is possible in a 
fraction of the Swift bursts. 
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1 INTRODUCTION 

Gamma-ray bursts (GRBs) emit extreme amounts of 7-rays 
on a short time scale; typically 10 50 — 10 54 erg are released in 
0.1 — 100 seconds. Only violent processes, such as a compact 
object merger or the collapse of a massive star, can explain 
these large energy releases, which have made GRBs observ- 
able out to high redshifts of z « 8 - 9 ilTanvir et al.l 120091 : 



ISalvaterra et al.ll2009l ; ICucchiara et al.ll201lh ~ 

There exists strong evidence, that the collapse of 
massive stars can pro duce long GRBs (Tgo > 2 s; 
iKouveliotou et alJlilm since spectroscopic features from 
supernovae (SNe) have been dete cted in optica l follow - 
up observat i ons o f GRBs (e.g. Galama et al.l (|l998h ; 
IHiorth et all (120031): IStarling et ah! (120111 ). Also see review 
by iHiorth fc Blooml ^Olll )). All these SNe are of type Ic 
with broad lines and no signs of Hydrogen or Helium. Be- 
sides these spectroscopic detections, evidence for SN Ic fea- 
tures is also found in light curves of some GRBs JPatat et al.l 
120011 : iLipkin et al.l2004l : ICobb et al.ll2010l : ICano et al.ll201ll ). 

One burst of particular interest wa s GRB 060218 
|Soderberg et al.ll200fi; iMazzali et al.ll2006l ). It had an asso- 
ciated supernova ( Pian et alj 20061 ). and its X-ray afterglow 
could best be described by a combination of synchrotron 
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emiss ion, which is usual for a fterglows, and b l ack bo dy emis- 
sion (jCampana et alj|2006l) . IWaxman et~aH (|2007l) showed 
that this black body emission could origin in a shock gen- 
erated by the breakout of a supernova through the sur- 
face of the GRBs progenitor star. Subsequently, thermal 
X-ray emission which may be described by a black body 
has been suggested in GRB 090618, GRB 10 0316D and 
GRB 1 01219B, which all have associated SNe dPage et all 
(|201ll ). IStarling et all (l201ll) IStarling et alj i|2012l ). respec- 
tively. See also lStarling et al.l (|2012l )). This supports the con- 
nection of the black body component with emission from a 
supernova. Deviations from a single power l aw in the earl y 
X-ray spectra in GRBs, was also found by iButlerl l|2007l ). 
who identified a soft emission component in 5 — 10 % of the 
bursts in the studied sample. 

In this series of papers, we search for more bursts with 
X-ray black body components, and derive the conditions 
under w hich such component s may be reliably recovered. In 
Paper I l|Starling et al.ll2012h black body components were 
identified in bursts with spectroscopic or photometric sig- 
natures in the optical. The aim of this paper is to perform 
a systematic search for more bursts with X-ray black body 
components in the Stmft sample, and to derive the conditions 
under which such components may be reliably recovered. In 
Section 2 the sample and model fitting are described, and 
in Section 3 we create simulated spectra to set constraints 
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on the detectability of black body components. In Section 4 
bursts are selected as candidates for having black body emis- 
sion, and Section 5 presents and discusses the final list of 
candidates. Section 6 extracts physical parameters, assum- 
ing that the excess emission is black body emission, and the 
fraction of GRBs with probable excess emission is examined. 

For the cosmological calculations we assume a ACDM- 
universe with h = 0.71, Q m = 0.27, and Q A = 0.73. All 
stated errors and error bars are 90% confident. In the plots 
rig is in units of 10 22 cm -2 unless stated otherwise. We will 
use the words thermal components and black body compo- 
nents interchangeably. 

2 FITTING MODELS WITH BLACK BODY 
COMPONENTS TO A SAMPLE OF GRBS 

2.1 A sample of Swift bursts 

The sample, in which we will search for X-ray black 
body compo nents, consists of the GRBs observed with the 
Swift XRT l|Burrows et al.l I2005T ), where redshifts are de- 
termined with optical spectroscopy, or simultaneous multi- 
band photometry, as o btained from the GROND-instrument 
(|Kriihler et al.ll2011al ). Furthermore we only select bursts, 
where Windowed Timing mode (WT mode) observations ex- 
ist, since this assures that data were taken shortly after the 
trigger. The sample consists of 190 bursts with the most 
recent burst from August 8th 2011. 

We use d the p ublicly available^ WT mode data, see 
I Evans et al.l ((2009). Note, that the sample also includes 
bursts not triggered by Swift. We make use of the time- 
averaged WT mode spectra, which have all been created in 
the same manner according to Evans et al. (2009) and us- 
ing the Swift software version 3.8 and the latest calibration 
data. 



2.2 Spectral modelling 

In order to identify black body emission for a burst, it 
is necessary to understand the interplay between all con- 
tributing components in the observed spectrum: the after- 
glow continuum, the black body emission and the absorp- 
tion of the source emission by gas in the line-of-sight within 
the host galaxy and the Galaxy. These properties can be 
parametrised by the seven parameters in Table[T] Note, that 
not all parameters are free: in our sample the redshift of 
each burst is known, and the Galactic colu mn density can, 
in mo st cases, be determined to within ±5% ( Kalberla et al. 
120051 ). 

For each burst we fitted four different models to the WT 
spectrum. These models are also shown in Table [T] Model 1 
(Ml) is an absorbed power law, whereas Model 2 (M2) also 
includes black body emission at the redshift of the bursts. 
Model 3 and 4 (M3 and M4) are similar to Model 1 and 
2, with the only difference, that riH.mt is fitted to the late- 
time Photon Counting mode spectrum, which we expect to 
be free from additional components and therefore a reliable 
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measure of any absorbing column, instead of from the early- 
time (WT) spectru m itself. The a tomic data used are So- 
lar abundances from Wilms et al. I <|2000l ) and cross-sections 
from IVerner et al.l (|l996T ). 

For all the WT spectra we group the spectra in 20 
counts per bin and perform a x 2_n t- For the late time ob- 
servations, used to determine nH.int in Model 3 and 4, un - 
grouped spectra are fitted with Cash statistics ICashl l|l979h . 

2.3 Selecting candidates with the F-test 

The F-test (e.g.; lLampton et al.lll976T ) is a statistical test, 
which gives the probability that an improvement in the re- 
duced x 2 > with the inclusion of extra parameters, is due 
to an improvement in the fitting model. To test whether a 
black body component is present in a spectrum, we therefore 
us e the F-test to c o mpar e Model 1 and 2 and Model 3 and 
4. IProtassov et al.l l|2002h showed, that the F-test is not a 
stringent test, but in the following section we will explicitly 
show, that the F-test in many cases is sufficient to recover 
black body components, when they are present. In Section [4] 
we will select the candidates. 



3 RECOVERING A BLACK BODY 

COMPONENT IN MOCK SPECTRA OF 
GRB 101219B 

In Section [3] it will be found, that one of the candidates 
for having a black body component is GRB 101219B with a 
temperature of 0.22 keV in the Model 2 fit (see Section|4]and 
Table [5}. This value is consistent wi th the value found in Pa- 
per I . With a redshift of 2 = 0.55 (|de Ugarte Postigo et al.l 
l201ll) GRB 101219B is one of the most distant GRBs wit h 
a spectroscopically detected supernova dSparre et al.ll201ll ). 
Now mock spectra will be used to derive the range of col- 
umn densities, redshifts and afterglow parameters for which 
a detection of the black body component is possible, and we 
will use this information to refine our candidate list. 

3.1 The role of column densities 

To reveal the role of the column densities, the F-test signifi- 
cances of the black body detections are calculated for several 
101219B-like mock spectra with nn-values varying from 0.0 
to 0.8 x 10 22 cm -2 . The other parameters used to create the 
mock spectra are identical to the recovered value in the fit 
of Model 2 from the real data. 

The left panel in Figure [T] shows how the F-test sig- 
nificance (calculated from Model 1 and 2) depends on the 
n_fj-value used to generate the mock spectra. If the F-test 
significance is lower than 10~ 4 we will say that Model 2 is 
favoured over Model 1 (i.e. the presence of a black body com- 
ponent is favoured). The central and the right panel show 
the tiH and kT recovered by the fits of the same mock spec- 
tra. When the uh used to generate the mock data is larger 
than 0.4 x 10 22 cm" 2 , the nu recovered by the fit is often too 
small, and in most cases a temperature with a large error is 
recovered. For nu < 0.4 x 10 22 cm -2 the F-test favours the 
existence of a black body component, and correct kT and 
n h -values are recovered within their errors. 

Similar mock spectra without black body components 
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Table 1. The parameters in the four fitting models (M1-M4). It is shown whether a parameter is free (if marked as free), fixed {fixed), 
not included (-) or fitted to the late time spectrum (late fit). 





Description 


Ml 


M2 


M3 


M4 


z 


Redshift of the GRB 


fixed 


fixed 


fixed 


fixed 


n H,gal 


Galactic column density 


fixed 


fixed 


fixed 


fixed 


"H.tat or n H 


Intrinsic column density in the GRB host galaxy 


free 


free 


late fit 


late fit 


r 


Photon index, defined as Fp oc E~ r (Fg is flux) 


free 


free 


free 


free 




Flux of the power law component of the afterglow in a 0.3 — 10 keV 


free 


free 


free 


free 




band 










kT 


Black body temperature (in the rest frame of the burst) 




free 




free 


Fbb 


Flux of the black body emission in a 0.3 — 10 keV band 




free 




free 




0.0 0.2 0.4 0.6 0.8 0.0 0.2 0.4 0.6 0.8 1.0 °' 1 0.0 0.2 0.4 0.6 0.8 

intrinsic nu (faked) intrinsic nj (faked) intrinsic uh (faked) 



Figure 1. Left: The F-test significance of the recovery of a 101219B-like black body component as function of nn (in units of 10 22 
cm -2 ). Central: The n^f-value recovered (in the Model 2 fit) as function of the correct njj- Right: The recovered temperature (in keV). 



were also made, and no significant detections of black body 
components were favoured by the F-test. 

3.1.1 Fixing nn to the late time spectrum 

For the fits of the real data the n_ff-values in Model 3 and 
4 are fixed to the fitted value from the late time spectrum. 
Now it will explicitly be shown, that this method might lead 
to significant black body detections for spectra with no real 
black body components, when uh is fixed to a too large 
value. 

First spectra with photon indices and power law nor- 
malizations identical to GRB 101219B are generated, but 
with nn varying from to 0.8 x 10 22 cm -2 . All the spectra 
are generated without black body components. In the fitting 
function nu is fixed to 0.5 x 10 22 cm -2 for all the spectra, 
just like if this value had come from the fit of a late time 
spectrum. In Figure [2] (left panel) it is shown that the F-test 
favours a black body component when nn < 0.35 x 10 22 
cm -2 . All the recovered black body components had tem- 
peratures between 0.08 keV and 0.2 keV. We conclude that, 
if uh is fixed to a value ~ 0.15 x 10 22 cm -2 larger than 
the real value, a spurious black body component will com- 
pensate for the flux lost due to the high column density 
absorption parameter. 

A side result of this analysis is that typical Galactic 
column densities can not lead to spurious detections of black 
body components. nH, ga i is typ ically of order ~ 0.05 x 10 22 
cm -2 with a 10 per cent error l|Kalberla et alJfeoOSl ). which 
is well below the critical value of 0.15 x 10 22 cm -2 found 
above. 



3.2 Redshift and afterglow parameters 

To find the redshifts for which it is possible to observe a 
GRB 101219B-like black body component, simulations like 
those in Section \'S. II are performed, with the only difference 
that z is varied instead of nn- Figure[2]( centra/ panel) shows, 
that black body components are recovered with F-test signif- 
icances < 10 -8 when z < 1. For 1 < z < 2 the significance of 
the detections are lowered by several orders of magnitudes. 

The role of the parameters describing the synchrotron 
emission from the afterglow will now be quantified. First we 
fix the photon index to the GRB 101219B- value and simulate 
spectra with power law normalizations covering the values 
observed in our sample. The resulting F-test significances of 
our fits are shown in Figure [2] (right panel). A low normal- 
ization constant is essential in order not to dilute the signal 
from the black body component with noise from the after- 
glow emission. We also fitted simulated spectra with a fixed 
normalisation and a varying photon index, see Figure [3] A 
trend in the complicated pattern is that GRB 101219B-like 
black bodies are most likely to be detected with low photon 
indices. 



4 SELECTING CANDIDATES 

Now we will start searching for candidates with black body 
components in our sample. The bursts with F-test probabil- 
ities lower than 0.005 in the Model 1-2 comparison or the 
Model 3-4 comparison were selected. For now, we are only 
interested in making an initial list of candidates, which will 
be refined later. We are therefore using a looser F-test sig- 
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Wrong fix of n# The role of redshift Power law normalisation 
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intrinsic n H (faked) z log 10 PLNORM 



Figure 2. Left: In the fits the n#-value is fixed to 0.5 X 10 22 cm -2 , which is different from the njj used the generate the spectra, shown 
as abscissa. Spurious detections appear, when uh 1s fixed to a to large value. Central: The F-test significance as function of z for a 
GRB 101219B-like burst. Right: The role of the power law normalization (in arbitrary units). 



nificance threshold here than in Section [3] (a limit of 0.005 
is now used instead of 10~ 4 ). The selected candidates are 
listed in Table 

4.1 Ruling out candidates with uncertain uh 

It is important to understand the contribution of ns in order 
to make a reliable detection of a black body component. We 
therefore removed all bursts with a 90% error on uh larger 
than 0.4 x 10 22 cm" 2 for Model 2 or larger than 0.8 x 10 22 
cm -2 for Model 4. A lower threshold was chosen for Model 
2, so candidates with a degenerate n //-value were excluded. 
In Table [2] a burst is marked with uncertain uh, if it is 
excluded due to this criterion. The analysis was done for 
larger values of the thresholds, but it did not lead to the 
detection of more candidates in the refined analysis, which 
will be presented in Section T4.3I and T4.4I 

4.2 Candidates at high and low redshift 

4-2.1 The high-redshift candidates (z > 2) 

For z > 2 GRB 101219B-like black body components are dif- 
ficult to detect (see Section 13.21 and Figure [2{- A number of 
candidates with z > 2 have, however, emerged in our analysis. 
They typically exhibit high temperatures (kT > 1 keV), which 
is much higher than most previously c laimed examples (e.g.; 
ICampana et alj2006l : IStarling et~aT]|201lh . They also have an un- 
certain nui which increases the risk of a spurious detection. 

At low redshift (e.g. z < 1), only a few candidates have 
temperatures larger than 1 keV, even though a detection of such 
a black body component is expected to be easier at lower redshift. 
The high-redshift candidates are therefore unlikely to represent 
real black bodies, and they will not be examined further in this 
work. 

4-2.2 The low-redshift candidates (z < 2) 

Three of the candidates, GRB 060218, GRB 090618 and 
GRB 101219B, already have identified black body components 
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Figure 3. The F-test significance as function of the photon index. 



in other studies, and they also have a spectroscopically confirmed 
SN in the optical. Among our candidates is also GRB 100418A, 
which also has a SN in the optical (de Ugarte Postigo et al. in 
prep.). Paper I studied the early X-ray spectrum of GRB 100418, 
and also find that a model with a powcrlaw and a black body 
model can give a good fit. However, they cannot rule out an ab- 
sorbed power law model given the large uncertainty on the in- 
trinsic absorbing column of this GRB host, which may be higher 
than the limit for detectability we find in Section l3.1l of 0.4 X 10 22 
cm -2 . Paper I also studied GRB 081007, and found an indication 
of a black body component, but no conclusive evidence. 

At first s ight it is worrying, that GRB 100316D, which 
IStarling etafl ^Ollr ) showed had an excess black body emission, 
is not among our candidates. An examination of this burst shows, 
that it is not selected due to a F-test value of 0.007 (for the 
Model 1-2 comparison), which is slightly above our threshold of 
0.005 defined in Section 12.31 It is therefore clear that our candi- 
date selection algorithm will not find all bursts with black body 
components (discussed further in Section l6.H I. 

Figurc|4]and[5]show how the F-test significances of the black 
body component detections in a selection of our low redshift 
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Table 2. The initial candidates for having a black body component. Candidates with z > 3 are omitted. The F-test significances are 
marked as bold if the burst is a candidate in the given model. If the recovered nn is too uncertain the burst is marked with "una n^f" 
(see Section l4.1l for details). 



GRB 


2 


F-stat (Ml-2) 


F-stat (M3-4) 


kT (M2) 


050724 


0.257 


2.1 x 10~ 7 


overflow 


0.9 ±0.1 


050820A 


2.6147 


1.0 x 10~ 7 


overflow 


4+3.6 


060124 


2.3 


0.041 


3.4 X 10~ 8 


0.6 ±0.1 


060202 


0.783 


4.2 X 10~ 6 


4.3 x 10~ 14 


U.3»_ . 05 


060218 


0.0331 


1.0 X 10- 174 


9.6 X 10" 110 


0.123 ±0.002 


060418 


1.49 


3.0 X lO" 6 


overflow 


0.61 ± 0.05 


060502A 


1.5026 


5.6 x 10~ 5 


4.7 X 10" 5 


0.35 ± 0.04 


060604 


2.68 


0.16 


2.1 x 10~ 9 


n i c+0.09 
u - lo -0.02 


060714 


2.7108 


0.27 


2.0 x 10~ 9 


24+ ' 06 
u ' z4 -0.11 


060904B 


0.7029 


0.0014 


1.0 


n 9Q+ ' 06 

u ' ZM -0.05 


061021 


0.3463 


2.9 X 10~ 7 


3.2 X 10 15 


u - ±z -0.02 


061110A 


0.7578 


6.5 X 10 11 


0.00019 


0.32 ±0.02 


061121 


1.3145 


1.5 X 10~ 8 


1.1 X 10~ 9 


0.50 ±0.06 


070318 


0.8397 


0.0075 


9.9 X 10- 11 


0.3 ±0.1 


070419A 


0.9705 


0.17 


1.2 x 10~ 7 


0.4 ±0.1 


070508 
070724A 


0.82 
0.457 


1.0 x 10" 5 
0.15 


0.074 
0.0014 


0fi6 +0 ' 002 

U.UDD_ 004 
u ' lo -0.04 


071031 


2.6918 


6.1 x 10~ 9 


8.4 x 10~ 67 


2 1 +0 ' 4 

-0.3 


071112C 
080210 


0.8227 
2.6419 


0.00011 

0.10 


0.00013 

0.0020 


40+ 06 

1 n +0 ' 3 

A " u -0.2 


080310 


2.4274 


0.052 


1.7 x 10" 22 


z ' u -0.7 


080319B 
080413A 


0.9382 
2.433 


1.4 x 10" 25 

0.27 


1.0 

2.0 x 10~ 7 


1 02+ ' 08 
A - uz -0.07 
1 4+0.04 
u - ly: -0.03 


080430 


0.767 


0.00020 


3.3 X 10" 5 


0.24 ± 0.06 


080603B 


2.6892 


0.042 


2.8 x 10~ 13 


_i_n « 

A -'-o.3 


080604 


1.4171 


0.0010 


3.0 X 10 6 


0.34+™° 


080605 
080721 


1.6403 
2.5914 


0.00011 

2.6 X 10~ 10 


0.84 
5.3 x 10~ 51 


0.087±°;°g> 

21 + 01 
u - zl -0.02 


080805 


1.5042 


0.38 


1.1 x 10~ 7 


06+ 05 


080928 


1.6919 


2.9 X 10" 5 


overflow 


i 7 +0.4 


081007 


0.5295 


0.0016 


0.00032 


20+ 0,07 
u ' zu -0.05 


081008 


1.967 


0.057 


6.7 X 10 33 


2.1 ±0.8 


081109 


0.9787 


0.047 


1.0 X 10~ 5 


1 o +0 - 3 


081203A 


2.05 


0.00073 


0.053 


0.5 ±0.08 


081222 
081230 


2.77 
2.03 


0.086 
0.12 


2.9 x 10~ 5 
0.00023 


lfi+0 05 

9 S+0.7 
z - 8 -1.0 


090418A 


1.608 


0.062 


8.3 x 10~ 5 


0.5 ±0.2 


090424 


0.544 


0.0041 


0.0042 


i >s+0- 2 


090618 


0.54 


0.00026 


4.1 X 10~ 6 


1.9 ±0.2 


090809 


2.737 


0.020 


4.3 x 10~ 7 


0.21 ± 0.03 


090812 


2.452 


1.5 x 10~ 5 


1.9 x 10~ 8 


i r+0.5 
A -°-0.2 


090814A 


0.696 


1.6 x 10- 10 


1.3 x 10~ 10 


0.28 ±0.03 


090926B 


1.24 


2.9 x 10~ 9 


1.8 x 10" 9 


0.53 ± 0.07 


091029 


2.752 


2.0 x 10~ 5 


0.0054 


n Qo+0.08 


100117A 


0.92 


0.0055 


1.7 x lO" 10 


1 2+°- 4 


100418A 


0.6235 


0.00074 


2.3 x 10 11 


21 + 03 
u ' zl -0.04 


100621A 


0.542 


0.00011 


2.7 x 10~ 5 


U-48_o.oe 


101219B 


0.55 


2.2 X 10 13 


3.9 x 10" 27 


0.22 ± 0.03 


110205A 


2.22 


0.021 


4.7 x 10~ 10 


3±1 


110715A 


0.82 


0.40 


6.6 X 10~ 6 


ln+ 006 


110731A 


2.83 


0.22 


0.0049 


12+ 007 
u - iz -0.04 


110808A 


1.348 


0.0017 


0.00068 


0.29 ± 0.07 



kT (M4) 


n H (M2) 


n H (M4, fixed) 


Note 


underflow 


0.33 ±0.06 


0.2 ±0.1 




underflow 


Q +20 
-0 


00+ 09 

U.UU_Q QQ 


unc. riu High-2 


0.46 ±0.07 


0.6 ±0.1 


0.8 ±0.2 


High-2 


0.34 ± 0.03 


1.6 ±0.1 


1.7±0.2 




0.186 ±0.002 


0.76 ±0.02 


0.40 ± 0.04 




00+ 03 


0.5 ±0.1 


2 +0 - 3 
u ' z -0.2 




0.31 ±0.02 


0.4 ±0.2 


0.5 ±0.2 


unc. riu (M2 only) 


A - D -0.2 


z - z -0.4 


1.1 ±0.3 


High-2 


1.9 ±0.3 


3.0 ±1.1 


i 1+0.6 
1 - i -0.5 


unc. nil High-2 


n no+ 0,05 
u - uu -o.oo 


n co+o.07 

u -°'-'-0.06 


0.4 ±0.1 




0.12 ±0.01 


0.08l°; I 


0.08 ±0.02 




overflow 


0.14 ±0.04 


3+ - 4 




0.45 ±0.03 


0.5 ±0.1 


0.61 ±0.08 




0.24 ± 0.04 


r, ^ + 0.2 

°- 6 -o.i 


0.8 ± 0.1 




0.18 ±0.04 


0.5 ±0.1 




unc. riu 


0-030l°;°oi 


„ _ „_i_n oq 
A ' A -0.07 


0.6 ±0.2 




0.8«; 2 


1.0 ±0.5 






1.44 ±0.06 


0.9 ±0.1 


o n+ 
u - u -o.o 

r. 1 +0.2 

O-i-o.i 


High-2 


„ o ~+0.03 
°- 37 -0,04 


n r.o+0.08 

°- 08 -0.07 




i n +0 ' 3 




i c+0.8 
A - £) -0.7 


unc. riu High-2 


1.3 ±0.1 


kq+0.07 


4+°- 3 


High-2 


underflow 


0.12 ± 0.01 


^ ^o+0 04 

- 08 ±o:o3 




1.0 ±0.1 


o q+0.9 
z - 8 -0.6 


U -°-0.5 


unc. riu High-2 


0.21 ±0.03 


q+0-3 


0.36 ± 0.06 




-i ^.+0 2 


0.8 ± 0.3 


0.2+^ 


unc. nil High-2 


0.28 ± 0.04 


- .+0 21 

°- 04 -o:o4 


n r»4-0 5 

0.2+^;^ 




04+ 09 


-] 9 + 0.2 

± - z -o.i 


0.6 ±0.3 




2.0 ±0.1 


2.2 ±0.3 


0.5 ±0.2 


unc. nil High-2 


0.27 ±0.05 


65+ 16 


1.2 ±0.5 


unc. nn 


underflow 


0.5 ±0.1 


0.3 ±0.1 




0.30 ± 0.03 


0.8±° 3 2 


0.5 ±0.1 


unc. nn (M2 only) 


1 in+0.08 


0.9 ±0.1 


0.3 ±0.3 




0.13 ±0.04 


0.5 ±0.1 


1.0 ±0.2 




overflow 


1.4 ±0.5 


0+2 
z -l 


unc. nil High-2 


1 O + 0.3 

1,o -0.2 


1 0+ ' 4 


0.4 ±0.2 


High-2 


1.5 ±0.4 


0-8±a3 


0.2 ±0.2 


High-2 


0.27 ±0.09 


o.5lS:S 


1.1 ±0.2 




1.3 ±0.2 


44+0- 03 
u '^-0.02 


0.42 ±0.06 




0-381^3 


0.26 ± 0.02 


0.19 ± 0.02 




2.1 ± 0.3 


6±1 


0.6 ±0.3 


High-2 


0.10 ±0.02 




1 2+°- 7 

A - z -0.6 


High-2 


0.26 ±0.01 


0.14 ±0.06 


2+°- 3 
u - z -0.2 




0.52 ±0.04 


i c+0.3 
1 - o -0.2 


1 5+ 12 


unc. nil 


overflow 


1 0+ ' 7 


0.5 ±0.2 


unc. nil High-2 


0.09 ±0.2 


0.2 ±0.2 


n 9+1-6 
z ' z -1.2 


unc. nil 


0.24 ±0.01 


0.2 ±0.1 


o i+°- 4 




29+ 03 


1.4 ±0.2 


2.0 ±0.2 




0.20 ±0.1 


06+ 06 


0.1 ±0.1 




1 6+ - 3 


54+ - 08 


0.4 ±0.1 


High-2 


15+ ' 04 


1.0 ±0.4 


1.3 ±0.3 




i 7 +0.6 


1.8 ±0.3 


1 l+°- 7 


unc. nil High-2 


n 99 +0.02 
u ' zz -0.03 


2+ ' 3 
u - z -0.2 


6+°- 4 


unc. nil 
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C t).0 0.2 0.4 0.6 0.8 1.0 

intrinsic njj (faked) 



Figure 4. Same as Figure [T] (left), but with three additional 
bursts. 




-16 1 1 1 1 1 

0.0 0.5 1.0 1.5 2.0 

z 



Figure 5. Same as Figure [2] (central) , but with three additional 
bursts. 



candidates depend on column density and redshift. In general 
the same trends found for GRB 101219B above applies for both 
GRB 090814A and GRB 061021. A black body component, simi- 
lar to the one present in GRB 060218, is expected to be detectable 
at high redshifts (e.g. z > 2). 



4.3 Approach I: Fitting with 5 free parameters 

The final step in our selection procedure is to do a time sliced 
analysis of the remaining candidates, i.e. those without remarks 
in Tablc[2] In this section we will focus on the candidates selected 
in the comparison between Model 1 and 2. In Section l4.4l we will 
fit and analyse the candidates selected with Model 3 and 4. 

For each of the candidates, where Model 2 is favoured by 
the F-test, we manually inspected the light curves, and selected 
time intervals, where the light curves exhibit a single power law 
decay. We especially tried to avoid time intervals with flares and 
plateaus. 

To each time sliced spectrum all the models were fit, see a 
summary of the fits in Table [3] Several candidates with no flare- 
free epochs are excluded from the table. The candidates with ex- 
cess emission are GRB 061021, GRB 061110A and GRB 090814A. 



They have low column densities (n^f < 0.3 X 10 22 cm -2 ) for both 
Model 2, where the n^j is a free parameter, and Model 4, where 
nii is fitted with the late time spectrum. In all the fits the pres- 
ence of a black body component is favoured by the F-test with 
significances from 0.0005 to 2.4 X 10 — 13 . The temperatures are 
between 0.1 and 0.3 keV, which is consistent with claimed black 
body components in other studies. 

In Table \3\ some bursts are also listed as candidates with 
some caveats. These are GRB 100621A (because it has a F-test 
significance of 0.00144, and a large njj), GRB 061121 (it has large 
kT and 2, and the F-test does not favour a black body component 
in Model 3 and 4) and GRB 050724 (it has an unexpectedly high 
kT) . We note that GRB 050724 has been clas sified as a short burst 
with extended emission ( Barthclmy ct al . 2005), so its origins may 
not lie in the collapse of a massive star as is the case for the rest 
of our sample which are all long GRBs. 



4.4 Approach II: Breaking Nh -degeneracies with 
the late time spectrum 

The three bursts, which were candidates for having black body 
components in the previous subsection, all have low column den- 
sities < 0.3 X 10 22 cm" 2 . This is expected from Figure [4] 
which showed, that it is hard to recover a black body component 
at large njj, when it is fitted as a free parameter (as in Model 1 
and 2). 

We will now do a new search for black body components, 
where we use the the late time spectrum to fit nn, i.e. use Model 
3 and 4. This method can break the degeneracies between uh 
and the other parameters, but it also has the disadvantage that a 
spurious black body component can be found if a too high value 
of rin is recovered in the late time spectrum (see Section 13.1.10 . 

We fitted Model 3 and 4 to the candidates. Table [4] sum- 
marizes the fits. Several of the candidates are not shown in the 
table, since they have complex and bumpy light curves, or too 
little data to extract a good spectrum. 

Several bursts have favourable F-test significances. To ex- 
plicitly check whether the fixing of nu might have caused spu- 
rious black body detections, we made additional fits, where uh 
was fixed to a range of different values within the 90% confi- 
dence intervals recovered in the late time fits. Figure|6]shows how 
the significances of the detections depend on the value to which 
nii is fixed. The detection is not robust for GRB 060202 and 
GRB 061121 because of degeneracies between kT and nu from 
the late time spectrum. The list of candidates from the Model 
3 to 4 comparison consists of GRB 061021, 061110A, 081109, 
090814A, 100621A and 110715A. 



5 THE FINAL CANDIDATES AND 
POTENTIAL CAVEATS 

We will refer to the bursts in Table|4]as our final list of candidates 
for having a black body component. In this section we will discuss 
potential caveats for each candidate. 



5.1 The Redshifts 

It is essential that the redshift, which is a fixed parameter in 
all our models, is reliably determined. Table [57TI summarizes how 
the redshift has been found for the bursts with excess emission. 
GRB 061021, GRB 090814A and GRB 110715 have redshifs found 
with only two absorption lines. It is however possible that these 
absorption lines comes from other galaxies than the host galaxy of 
the bursts. The remaining bursts have more reliably determined 
redshifts with several absorption and/or emission lines. 
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Table 3. The time sliced fits from Section 14,31 where bursts are selected with Model 1 and 2. Three reliably detected black body 
components are found. 



Candidates with probable excess emission: 



GRB 


z 


Time 


F-test Ml 


-2 


F-test M3-4 


kT (M2) 


kT (M4) 


n H (M2) 


n H (M4,fixed) 


061021 


0.3663 


87 - 170 


3.1 x 10~ 


7 


2.4 x 10" 13 


13 +(J.U3 


0.13 ±0.01 


n ns+ u - i;:! 


0.08 ±0.02 


061110A 


0.7578 


180 - 240 


2.5 x 10" 


6 


0.00050 


0.25 ±0.01 


0.19 ±0.01 


o oo +0 - 07 
u - uu -o.oo 




090814A 


0.696 


166 - 265 


4.0 x 10" 


7 


2.2 x 10" 9 


0.32 ±0.03 


0.30 ±0.01 


0.12 ±0.06 


2+ ' 3 
u - z -0.2 


0908 14A 


0.696 


265 - 390 


1.1 x 10" 


(i 


2.1 x 10" 8 


n iq+0-05 
u - iJ -0.04 


0.20 ±0.02 


2 +0 - 2 
u - z -o.i 


2+ - 3 
u - z -0.2 


Candidates with 


some caveats: 














GRB 


z 


Time 


F-test Ml 


-2 


F-test M3-4 


kT (M2) 


kT (M4) 


n H (M2) 


n H (M4,fixed) 



050724 0. 
061121 1. 
100621A 0. 
100621A 0. 
Candidates 
GRB z 



257 
3145 
542 
542 

where 



100 
125 



190 1. 
215 1. 
80 - 120 
190 - 230 
a black body 
Time F 



3 x lO -0 
1 x 10~ 5 

0.022 

0.0014 
component 
-test Ml-2 



0.9 ±0.1 
0.44j 
0.69"^ 



060202 


0.783 


148 


- 303 


0.0023 


0.12 


060202 


0.783 


745 - 


- 1000 


0.28 


0.0058 


060418 


1.49 


240 


- 450 


0.29 


0.20 


070508 


0.82 


200 


- 500 


0.19 


0.17 


070508 


0.82 


500 - 


- 1000 


0.055 


0.33 


071112C 


0.8227 


89 - 


- 179 


0.022 


undef 


080319B 


0.9382 


200 


- 500 


0.0077 


1.0 


080430 












080604 


1.4171 


125 


- 200 


0.14 


0.059 


080605 


1.6403 


300 


- 700 


0.054 


1.0 


080928 


1.6919 


280 


- 320 


0.0014 


undef 


090424 


0.544 


700 - 


- 1000 


0.62 


0.017 



1 

0.86 
6.0 x 10"' 
is not favoured: 

F-test M3-4 kT (M2) 



+0.05 
0.06 
+0.06 
0.07 

0.33 ±0.06 



0.36 ±0.09 
O.Oj^J 0.6 ±0.1 

0.05 ±0.05 1.4j 
0.25 ±0.03 1.3"* 



„+0.6 
-0.3 
i O+0.4 
1 -° -0.3 



0.2 ±0.1 
0.61 ±0.07 
2.0 ±0.2 
2.0 ±0.2 



042 TrnT7- 

u - 4z -0.06 
overflow 

0.6 ±0.2 
+0.09 

0.02 
+0.8 
0.2 
+0.08 
0.09 
+0.4 
0.2 



kT (M4) 



045+D.U03 
u.u<io_ 002 

25+ - 08 
u - zo -0.06 

6 +0 - 2 
0.031 ±0.002 



n H (M2) n H (M4,fixed) 



0.06 
1.0 

0.40 
1.0 



0.21] 



,+0.15 

-0.06 

n ns4+ ' 004 
O.U84_ Q 00g 

2.3+ 09 



0.030 
0.00 
0.00 

0.27 



+0.014 
0.004 
+0.04 
0.00 
+0.03 
0.00 



1.6 ±0.1 
1.6 ±0.2 
0.2 ±0.2 

1 o+ - 3 

0.8 ±0.2 
0.2 ±0.1 
0.13 ±0.03 



1.7 ±0.2 

1.7 ±0.2 

2 +0 - 3 
u - z -0.2 

0.6 ±0.2 
0.6 ±0.2 
1+ ' 2 
0.08 ±0.03 



+0.06 

0.07 



003+ ' 094 



0.3 1 
1.3j 



1.1 



0.6 
+0.8 
0.4 



1+0.1 
-0.0 



o.o: 

0.8 ±0.2 



3+0.4 
-0.3 
3+0.2 
-0.1 

0.7 ±0.3 
0.52 ±0.07 



2+ ' 5 
u - z -0.2 

0.6 ±0.2 
0.3 ±0.1 
0.42 ± 0.05 



Table 4. The time sliced fits from Section 14.41 where bursts are selected with Model 3 and 4. Three reliably detected black body 
components are found. 

Candidates with probable excess emission: 



GRB 


z 


Time 


F-test M3 


-4 


kT (M4) 


n H (M4, fixed) 


061021 


0.3463 


87 - 170 


2.4 x 10" 


111 


0.13 ± 0.01 


0.08 ±0.02 


061110A 


0.7578 


180 - 240 


0.00050 




0.19 ± 0.01 


n+0.3 
U.O_ 2 


081109 


0.9787 


90 - 200 


6.3 x 10" 


■7 


0.16 ±0.03 


1 o +0 - 2 


0908 14A 


0.696 


166 - 265 


2.2 x 10" 


■9 


0.30 ±0.01 


0.3 ±0.2 


090814A 


0.696 


265 - 390 


2.1 x 10" 


■8 


0.20 ±0.02 


0.2 ±0.2 


100621A 


0.542 


80 - 120 


0.86 




0.05 ± 0.05 


2.0 ±0.2 


100621A 


0.542 


190 - 230 


6.0 x 10" 


-7 


0.25 ± 0.03 


2.0 ±0.2 


110715A 


0.82 


97 - 501 


6.6 x 10" 


-6 


0.15 ± 0.03 


1.3 ±0.2 


Candidates with 


some caveats: 








GRB 


z 


Time 


F-test M3- 


-4 


kT (M4) 


n H (M4, fixed) 


060202 


0.783 


148 - 303 


0.12 




n 045+ ,u)U4 
u.u+o_ 002 


1.7 ±0.2 


060202 


0.783 


745 - 1000 


0.0058 




25+ - 07 
u - zo -0.06 


1.7 ±0.2 


061121 


1.3145 


125 - 215 


7.5 x 10- 


6 


0.43 ±0.03 


61 + 07 
u - Di -0.06 



Candidates where a black body component is not favoured: 

GRB 2 Time F-test M3-4 kT (M4) n H (M4,fixcd) 



071112C 0.8227 89 - 179 undef 0.00+^ 

080604 1.4171 125 - 200 0.059 0.3 ±0.1 

090424 0.544 700 - 1000 0.017 0.8 ± 0.2 



1 +u " ! 
u - 1 -o.i 

2+ - 5 
u - z -0.2 

0.42 ±0.05 



5.2 Light curves and spectral evolution 

Figure [7] shows the parts of the light curves which were used in 
the time sliced analysis of the final candidates. We tried to avoid 
achromatic and X-ray flares and effects from flattening of the light 
curves, when the time intervals were selected, but we note that 
evolution of the spectral parameters, the power law index and/or 
the black body temperature and normalisation, may be occurring 



during any of our spectra. When analysing time sliced spectra 
this risk is lowered, and we examined the hardness ratios in the 
XRT Repository in order to flag incidences of strong spectral 
evolution. GRBs 061021, 061110A, 090814A and 110715A have 
approximately constant hardness ratios in the intervals where 
time sliced spectra were analysed. GRB 081109 shows a slowly 
hardening spectrum, while GRB 100621A becomes significantly 
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Table 5. A schematic summary of how the redshift has been determined for each of the final candidates for having black body components. 



GRB 



Redshift Absorption lines 



Emission lines 



References 



061021 

061110A 

081109 

090814A 

100621 



0.3463 

0.7578 

0.9787 

0.696 

0.542 



110715A 0.82 



Mgn2796, Mgn2803 
Mgn2796, Mgn2803 

Mgn and Call 



Call, Cal 



H/3, [Om] doublet 
GROND photometry 

[On] 3727, H^, [Om] doublet 



Fynbo ct al 
Fynbo ct al 



Kruhlei ct al 



2009) 



2009) 



Jakobsson ct al. (2009 
Milvang- Jensen ct al. 



Kruhter et al 



Piranomonte et al. nBoilT) 



(2010): 



Fits inside the confidence interval of n# 




— *— 081109 



0.0 0.5 1.0 1.5 2.0 

n H 

Figure 6. In Model 3 and 4 nu is fixed to a value determined 
from late time observations. This figure examines how the F-test 
significances of the candidates from Table[4]depend on n^f-values 
with the 90% confidence intervals found in the fits of the late time 
spectra. 



softer during the first time slice but shows no spectral evolu- 
tion during our second time slice. Therefore the parameters we 
derive in Table [531 for these latter two GRBs should be consid- 
ered more uncertain than the error bars allow. Evolution of the 
black body component has been dem onstrated for GR Bs 060218 
(e.g. ICampana . et alj|2006h. 090618 llPage et alj|20ljl . 100316D 
dStarling et alj|201ll ; lOlivares E. et al.ll2012T) and 101219B (Pa- 
per I), where the black body cools and expands with time, but 
typically this evolution is slow compared with the exposure times 
covered by our X-ray spectra. 



5.3 Contribution from the prompt emission 

It is important that the prompt emission, which has a spectrum 
peaking at the energy, -E pea k> is n °t entering our X-ray spectrum 
for the Swift XRT (0.3 — 10 keV), since this potentially could give 
features resembling black body components. Here is a summary 
of the measurement of -E pea k f° r each burst: 

061021 Speak was meas ured to lie above 5 40 keV up to 8 sec- 
onds after the trigger ( Golcnctskii ct al . 2006). This is well above 
XRT-energies. 

061110A No spectral break or cutoff was seen with Swift /BAT, 
so no -Ep Ca k information is present. The role of the prompt emis- 
sion remains unknown for this burst. 

090814A No spectral break or cutoff was seen by Swift /BAT. 
The role of the prompt emission remains unknown for this burst. 



081109 Fermi/G BM measured £ poak = 240 ± 60 keV 
(von Kienlidl200Sl') after 26 seconds. Our spectrum is well below 
that energy. 

100621A Konus-Wind fits t heir spectrum covering — 74 sec- 
onds with _E peak = 95±J§ keV JGolenetskii et al.l2010T) . Our spec- 
trum starts at 80 s, so probably the peak energy was still well 
above the XRT band. 

110715A 18 second s after trigger Konu s-wind measured 
S pcak = 120+ j? keV jGolenetskii et ; alJl2uTll ). and BAT found 
120 ± 21 keV ilJkwatta et al.ll201lf) . It is unlikely, that this af- 
fects our XRT spectrum, which starts at 98 seconds. 

5.4 Is it a black body component? 

The six bursts in our final candidate list (GRB 061021, 
GRB 061110A, GRB 081109, GRB 090814A, GRB 100621A and 
GRB 110715A) are clearly special, since models with a black body 
component included give better fits than absorbed power laws. 
We do, however, not have conclusive evidence that the actual 
emission mechanism is a black body component. Other possible 
explanations could be a model with multiple absorption compo- 
nents at different redshifts, the prompt emission peak could be 
moving through the band, there could be strong spectral evolu- 
tion, or the determined redshift could be wrong. We have a total 
of 190 bursts in our sample, so some of them will likely be af- 
fected by such features. Due to these caveats we will not report 
discovery of black body emission, but instead we will report the 
finding of bursts with possible black body emission. 



6 PROPERTIES OF THE BLACK BODY 
COMPONENTS 

In this section, we will assume, that the actual mechanism be- 
hind the excess emission in the spectra of the six bursts is a sin- 
gle temperature black body component. Table 15.31 shows radius 
(assuming spherical symmetry) and temperature for these bursts 
together with values reported in other studies. For our final can- 
didates the radii are calculated from the fitted temperature and 
luminosity. In the calculation we assume a spherically symmetric 
black body, in thermal equilibrium, which can be described by 
the Stcfan-Boltzmann law, 

L(erg s" 1 ) = 1.105 X 10 29 X r 2 (m) X T 4 (keV) (1) 

Also note that the luminosities and radii in the table are lower 
limits, since only the photons in a 0.3 — 10 keV are included in 
the luminosity calculation. Due to these caveats, the values in the 
table should be seen as rough estimates only. 

The bursts have temperatures in the range, 0.1 — 0.3 keV, 
which is consistent with the previously proposed black body com- 
ponents in the time-averaged WT spectra of GRB 060218, 090618, 
100316D and 101219B. They also have luminosities consistent 
with these previously studied bursts, whereas the radii are sligthly 
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Figure 7. The light curves for each of the bursts with probable excess emission. The curve is red (and solid) in the intervals selected for 
a time sliced analysis. 



Table 6. Radii, temperatures and luminosities of the black body components of the final candidates. Candidates from this work and 
other works are shown. If a value is bold, it has been calculated using Stefan-Boltzmann's law (see Eq. Jl}), from the two non-bold 
values for a given candidate. If a £; so -value is underlined , it has been calculated from the fluence in the 15-150 keV band and it should 
only be seen as a lower limit. 



GRB 


Radius (m) 


kT (keV) 


L (10 4Y erg s" 1 ) 


ref 


2 


Time (s) 


T 90 (s) 


Eiso (erg) 


061021 


6.6 X 10 1U 


0.13 


1.4 


this work 


0.3663 


87 - 170 


46 


4.6 x 10 51 


061110A 


1.4 X 10 11 


0.19 


25.6 


this work 


0.7578 


180 - 240 


41 


2.9 x 10 51 


081109 


2.0 x 10 11 


0.16 


27.1 


this work 


0.9787 


90 - 200 


190 


1.8 x 10 52 


090814A 


3.1 x 10 10 


0.30 


8.4 


this work 


0.696 


166 - 265 


80 


2.8 x 10 51 


090814A 


4.7 x 10 10 


0.20 


3.7 


this work 


0.696 


165 - 390 


80 


2.8 x 10 51 


100621A 


8.7 x 10 10 


0.25 


32.4 


this work 


0.542 


190 - 230 


64 


2.8 x 10 52 


110715A 


2.7 X 10 11 


0.15 


43.0 


this work 


0.82 


97 - 501 


13 


4.1 x 10 52 


060218 


1.0 x 10 1U 


0.20 


0.2 


Campana et al. (2006) 


0.033 


200 


2100 


6.2 x 10 49 


090618 


3.0 x 10 10 


1.00 


1000.0 


Page et al. (2011) 


0.54 


150 


113 


2.5 x 10 53 


090618 


1.0 x 10 11 


0.20 


100.0 


Page et al. (2011) 


0.54 


250 


113 


2.5 x 10 53 


100316D 


3 x 10 10 


0.14 


0.3 


Starling et al. (2011) 


0.0591 


250 


> 1300 


3.9 x 10 49 


101219B 


2.4 X 10 10 


0.20 


1.0 


Paper I 


0.5519 


250 


34 


4.2 x 10 51 



larger for GRB 061110A, 081109 and 110715A than the previously 
claimed examples. 

In t he case of GRB 060218 sev eral studi es (Ghiscllin i et al] 

l2007al lbl; IChevalier fc Franssonll200g|; IlH |2007|) disfavour the sce- 
nario that the excess emission is due to a black body component, 
since they can not explain the large black body luminosities. All 
our bursts with possible black body emission have a luminosity 
larger than GRB 060218 (Table [B~3l . so they might suffer from a 
similar problem. 

All the bursts have durations (Tgo) between 13 and 190 sec- 
onds, which is typical for long bursts. The redshifts are in the 
range, z = 0.37 — 0.98, so they arc slightly more distant than the 
cases reported in other studies (z = 0.03 — 0.55). 



6.1 The fraction of bursts with black body 
components 

In the Model 1 to 2 comparison we found three new candi- 
dates (GRB 061021, 061110A, 090814A) and we re-discovcrcd 
GRB 060218, 090618, 101219B. So we found a black body com- 
ponent in 6 out of the 116 bursts (i.e. 5%) with successful fits. 

We will now shed light on the fraction of Swift bursts for 
which it is possible to recover a black body component like the one 
in GRB 101219B. In Section|3]it was found that such a black body 
component can be detected (with the F-test comparison of Model 
1 and 2) if z < 1.5, tih < 0.4x 10 22 cm -2 and the logarithm of the 
powerlaw normalisation is below —0.5 (see Figure See Figure [l] 
[2] and [3)1 ■ Figure [8] shows nu versus z for all the bursts in our 
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Figure 8. njj versus z for all the successfully fitted bursts. If a 
black body component, like the one present in GRB 101219B, is 
detectable in a GRB it is marked with a black circle. In 15 out of 
102 bursts it is possible to recover such a black body component. 



sample with successful fits. It is marked whether a bursts passes 
these detection criteria. We see that in 15 %, 15 out of 102 bursts, 
a detection will be possible. This calculation is, of course, only a 
rough estimate. The X-ray black body components need not all 
be like the one in GRB 101219B, and the limits we have used 
on z, nu, r and power law normalisation are expected to be 
degenerate with each other. This calculation is, however, sufficient 
to establish that black body components only are detectable in 
a fraction of the bursts, even if such a black body component is 
present in all GRBs. 

None of the new candidates found in this work have a opti- 
cally confirmed SN. Note, however, that we excluded bursts which 
already arc discussed in other works, including Paper I where sev- 
eral bursts with associated SNe are analysed. 



7 CONCLUSION 

We have first examined under which conditions a black body com- 
ponent, like the one proposed to be in GRB 101219B, can be re- 
covered. At high rcdshift and in environments with intermediate 
or high column densities a detection of such a component will not 
be possible. We also find that it will be hard to recover a black 
body component, when a bright afterglow emission is present. 
We show that detection of a black body component only will be 
possible in a small fraction of all GRBs in our sample. 

We have searched for black body components in all the Swift 
bursts with known rcdshift and created a list of bursts with pos- 
sible black body components (GRB 061021, 061110A, 081109, 
090814A, 100621A and 110715A). They have temperatures, radii 
and luminosities similar to those found in previous studies of black 
body components in GRBs. 
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